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Abstract 

We present an algebraic approach to one-loop tensor integral reduction. The integrals 
are presented in terms of scalar one- to four-point functions. The reduction is worked 
out explicitly until five-point functions of rank five. The numerical C++ package PJFry 
evaluates tensor coefficients in terms of a basis of scalar integrals, which is provided by 
an external library, e.g. QCDLoop. We shortly describe installation and use of PJFry. 
Examples for numerical results are shown, including a special treatment for small or 

H vanishing inverse four-point Gram determinants. An extremely efficient application of 

the formalism is the immediate evaluation of complete contractions of the tensor integrals 
• • with external momenta. This leads to the problem of evaluating sums over products of 

. signed minors with scalar products of chords. Chords are differences of external momenta. 

^ These sums may be evaluated analytically in a systematic way. The final expressions for 

the numerical evaluation are then compact combinations of the contributing basic scalar 
functions. 



1 Tensor reductions 

The reduction of tensorial Feynman integrals to scalar Feynman integrals is an old technique. 
A systematic approach has been developed for one- to four-point functions in pp. The simple 
case of a self-energy vector 
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is evaluated in terms of scalar functions A , B as follows: 
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This works fine for many situations, but for higher-point functions and for specific kinematical 
situations other approaches are more useful. We advocate here Davydychev's approach [2] 
and represent n-point tensors by n-point scalars in higher dimensions and with higher indices 
(powers of the propagators Di = l/[(k — qi) 2 — m 2 }). An example is: 
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1 + Sij) (1 + Sik + 5jk)- The ig is defined here in D = d — 2e + 2L dimensions and 
propagators of lines i, j, k have a power raised by one unit. In a next step Tarasov's dimensional 
recurrences pH H] may be applied in order to diminish the dimensions of the scalar integrals]^] 
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where the operators i 1 * 1 , j^, act by shifting the indices z/j, z/j, Z/& by ±1. 

After dedicated simplifications [5j [6] , the tensor of rank R = 3 may be written as follows: 
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At this stage, we have reached: 
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4 The explanation of notations may be found in any of the articles |3J [U [SJ.For Gram determinants and 
signed minors see[A| 
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All scalar 5-point integrals are eliminated. 

No inverse Gram determinants Q5 have appeared. 

Scalar 4-point integrals in higher dimensions and/or with higher indices are present yet, 



e.g. I. 



[d+] 2 ,s 
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Inverse Cayley determinants (°) have appeared. 



The integrals 1^3 ' s arise from 5-point integrals by shrinking line s. 

By further, nontrivial manipulations [6] we have avoided shifted indices and isolated the 
inverse Gram determinants 04 in the higher-dimensional scalar 4-point functions; see e.g. in 
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These equations contain yet the generic 4-point and (partly indexed) 3-point functions in higher 
dimensions, ; /j^- , etc. Several strategies are now possible: 

• Evaluate them analytically in d + 21 — 2e dimensions - if you may do that. 

• Evaluate them numerically in d + 21 — 2e dimensions. 

• Reduce them further by recurrences — > apply ^ - but buy some towers of 1/04 at this 
stage. 

• Make a small Gram determinant expansion — » apply d6|) the other way round 



The last two items are done here. 

For non-small 4-point Gram determinants, the direct, iterative use of @ yields 
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and we are done. This works fine if 04 is not small [and also the Qj- 
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For small () 4 , one may apply the recurrence relation the other way round and get: 
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This works in fact fine, especially when improved by Pade approximations. For explicit ex- 
amples see [6] and sec. [2] on the P JFry package here. The approach has been independently 
implemented and used by Reina et al. [7]. 

2 PJFry 

The goal of the C++ package PJFry is a stable and fast open-source implementation of tensor 
reduction, suitable for any physically relevant kinematics]^] It performs the reduction of 5-point 
1-loop tensor integrals up to rank 5. The 4- and 3-point tensor integrals are obtained as a 
by-product. Main features are: 

• Any combination of (real) internal or external masses 

• Automatic selection of optimal formula for each coefficient 

• Leading 05 are eliminated in the reduction 

• Small 04 are avoided using asymptotic expansions where appropriate 

• Cache system for tensor coefficients and signed minors 

• Interfaces for C, C++, FORTRAN and Mathematica 

• Uses QCDLoop [HI E] or OneLOop pi)] for 4-dim scalar integrals 



Available from the project webpage https://github.com/Vayu/PJFry/ [HIE] 



The installation of PJFry may be performed following the instructions given at the project 
webpage. 

The project subdirectories are 
./src - the library source code 
./mlink - the MathLink interface 

./examples - the FORTRAN examples of library use, built with make check 

A build on Unix/Linux and similar systems is done in a standard way by sequential per- 
forming ./configure, make, make install. See the INSTALL file for a detailed description of the 
./configure options. 



5 The presentation relies partly on: V. Yundin, One loop tensor reduction program PJFRY, talk at meeting 
of SFB/TR9, 15 Nov. 2011, Aachen, Germany. 



The PJFry is used as one option of the GoSam package [12"] . 

The functions for tensor coeffcients for up to rank-5 pentagon integrals are declared in the 
Mathematica interface: 

In: = Names ["PJFry '*"] Names [PJFry] 

{AOvO, BOvO, BOvl, B0v2, COvO, COvl, C0v2, C0v3, \ 
DOvO, DOvl, D0v2, D0v3, D0v4, EOvO, EOvl, E0v2, \ 
E0v3, E0v4, E0v5, GetMu2, SetMu2} 

Out= {AOvO, BOvO, BOvl, B0v2, COvO, COvl, C0v2, C0v3, \ 
DOvO, DOvl, D0v2, D0v3, D0v4, EOvO, EOvl, E0v2, \ 
E0v3, E0v4, E0v5, GetMu2, SetMu2} 

The syntax is very close to that of e.g. LoopTools/FF: 

E0v3[i, j ,k,pls,p2s,p3s,p4s,p5s,sl2,s23,s34,s45,sl5,mls,m2s,m3s,m4s,m5s,ep=0] 

where: 

i, j, k are indices of the tensor coefficient (0<i<j<k<n), 

pls,p2s, ... are squared external masses p?, 

sl2, s23, ... are Mandelstam invariants (pi +pj) 2 , 

mls,m2s, ... are squared internal masses mf, 

ep = 0, — 1, — 2 selects the coefficient of the e-expansion. 

The average evaluation time per phase-space point on a 2 GHz Core 2 laptop for the eval- 
uation of all 81 rank 5 tensor form-factors amounts to 2 ms. 

Numerical examples are shown in figures [T] and [2] for a rank R = 4 tensor coefficient in a 
region, where one of the 4-point sub-Gram determinants becomes small 

£3333(0, 0, -6x 10 4 (x + 1), 0, 0, 10 4 , -3.5x 10 4 , 2x 10 4 , -4x 10 4 , 1.5x 10 4 , 0, 6550, 0, 0, 8315) 

When x — 0, the 4-point Gram determinant vanishes; see [A] for details. 



3 Contractions of tensor integrals with external momenta 



The contents of this section is worked out in detail in [13] . A starting point are representations 
of tensors in terms of basis functions I®, which are independent of the tensor indices {i, j, k, ...}. 
Then, one may use contractions with external momenta in order to perform all the sums over 
{i, j, k, ...}. This will evidently lead to a significant simplification and shortening of calculations. 
To demonstrate the idea, let us look at the most involved 5-point tensor studied: 
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In the approach avoiding inverse Gram determinants I/O 5, we have 
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Figure 2: Relative accuracy of S3333 near the region of small Gram determinant. At x ~ 0.0015 
PJFry switched to the asymtotic expansion. 



and simpler expressions. Now, in a next step, one may avoid the appearance of inverse sub- 
Gram determinants 04 and eliminate the higher indices of propagators {i,j, k, . . .}. Then, the 
complete dependence on the indices % of the tensor coefficients is contained in the pre-factors 
with signed minors: the tensor indices decouple from the integrals. 
As an example, we reproduce the 4-point part of I^^lf- 
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2 4 ( ' 



In (17), the 3-point terms are not reproduced, and one has to understand the 4-point integrals 
to carry the corresponding index s and the signed minors are (° k ) — > (j^) 5 etc. 

A chord is the momentum shift of an internal line due to external momenta, Di = (k — 
qi) 2 — mf + ie, and = {p\ +P2+ +Pi), with q n = 0. The tensor 5-point integral of rank R = 1 
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yields, when contracted with a chord, 
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and we get immediately the desired compact result for the contraction of chords (or 
momenta) and tensor integrals: 
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Similarly, the tensor 5-point integral of rank R = 2 may be treated: 
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It has the following tensor coefficients free of l/()s: 
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Equation (21) yields for the contractions with chords: 
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and finally (24) simply reads 
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Figure 3: Momenta definitions. 
This can be extended also to higher ranks. We need at most double sums, e.g.: 
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Many of the sums over signed minors, weighted with scalar products of chords are given in [13] , 
and an almost complete list may be obtained on request from the authors. 



A Determinants 

The Gram determinant Gn of the diagram shown in [3] is 

G N = \2qiqj\,i,j = 1,...N. 



(28) 



We are also using the modified Cayley determinant ()n of a diagram with N internal lines. For 
the choice qN = 0, both determinants are related: 
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where the matrix elements are defined by 

Yij = -(qi-qj) 2 + m 2 i + 
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Evidently, the Gram determinant Gn does not depend on the masses, and so doesn't ()#. 

Signed minors of Qjy are constructed by deleting m rows and m columns from Qn, and 
multiplying it with a sign factor: 
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(31) 



where sgn u\ and sgnr fe i are the signs of permutations that sort the deleted rows j± - ■ ■ j m and 
columns k\ ■ ■ ■ k m into ascending order. 
Example: 
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To be definite, take fig.[4]as a starting point. The example is taken from [Hj. The corresponding 
4-point tensor integrals are, in LoopTools/FF [151 HS1 [9] notation: 



D0i(id, 0, 0, s 9u , ted, t St M, s^p u , 0, Mf , 0, 0). 
The Gram determinant is 
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and it vanishes if: 
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In terms of a dimensionless scaling parameter x, 

ted (1 "I - 3^)^erf,crit) 

the Gram determinant becomes: 
The modified Cayley determinant 
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has another dimension. From (|14|) we see that a small Gram determinant expansion will be 

(39) 



useful when the following dimensionless parameter becomes small: 
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where S is a typical scale of the process, e.g. we will choose S = s 
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Figure 4: A six-point topology (a) leading to four-point functions (b) with realistically vanishing 
Gram determinants. 
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